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We suggest a closed form expression for the path integral of quantum transition amplitudes. We 
introduce a quantum action with renormalized parameters. We present numerical results for the 
V ~ X*' potential. The renormalized action is relevant for quantum chaos and quantum instantons. 
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;_! 1. Introduction 

O^The path integral has become a standard method to 
•^ quantize classical theories. The class of systems for which 
-^ the path integral can be computed analytically is ridicu- 
,__|lously small. Examples are quadratic Lagrangians, e.g., 
for free motion and the harmonic oscillator [n|. Let us 
Cn consider the Q.M. transition amplitude from Xin, tin to 
_^^/*> tfr given by 
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J dm^x^ — V[x) denotes the classical action. 



5^where S — , ^^ „ . 

■!i^In some cases this path integral can be expressed as a 

T^ sum over classical paths only 
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.5^ where S'[a;c!] is the classical action evaluated along the 

^^ classical trajectory from Xin, tin to Xfi, tfi. There may be 

5-H infinitely many classical paths for a given pair of bound- 

. ..ary points (xi„,ii„), {xfi,tfi) [|l],^. Different paths then 

correspond to different values of the action. The factor Z 

represents some (time-dependent) normalisation factor, 

enforced by the unitarity of the amplitude. Eq.(0) holds, 

e.g., for quadratic Lagrangians, where the sum runs over 

a single path where the action is minimal. In the case of 

the harmonic oscillator V{x) = ^io^x^ , one has 



group manifold of SU{2) 0]. Eq.(S) is a well known but 
remarkable result. It is possible that there are other sys- 
tems, where the path integral is given by a sum over clas- 
sical paths. However, the result is not true in general. A 
simple counter example is given by Nelson W. 

In the following we will explore if the validity of Eq. (||) 
can be extended to a wider class of quantum systems, if 
we allow the classical action to be replaced by a renor- 
malized action. The motivation comes from perturbative 
renormalisation in Q.F.T. Recall that a theory is called 
renormalizable, if a fixed finite number of counter terms 
of the action allows to make the Q.F.T. finite to any given 
order. The renormalized action is different from the bare 
action, due to quantum loop corrections, which exist for 
interacting theories. The action and also Green's func- 
tions have the same structure in its bare and renormal- 
ized form. Quantum mechanics can be viewed as Q.F.T. 
in + 1 dimension. Here we suggest for Q.M. the exis- 
tence of a renormalized/quantum action, which describes 
the transition amplitude (Green's function), having the 
same structure but parameters different from the classi- 
cal action. 



Conjecture: For a given classical action S = J dt^i 



V{x) there is a renormalized action S = J dt-. 
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which allows to express the transition amplitude by 
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Another example where the sum over classical paths is 
exact is the path integral of the quantum mechanical top, 
which mathematically corresponds to free motion on the 



Here id denotes the classical path corresponding to the 
action S, such that the action S{xci) is minimal (we 
exclude the occurrence of conjugate points or caustics). 
Z denotes the normalisation factor corresponding to S. 
Eq.(^ is valid with the same action S for all sets of 
boundary positions Xfi, xtn for a given time interval 
T ~ tfi — tin ■ The parameters of the renormalized action 
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depend on the time T. The renornialized action con- 
verges to a non-trivial limit when T —^ oo. Any depen- 
dence on xfi, Xin enters via the trajectory Xd- Z depends 
on the action parameters and T, but not on Xfi , Xin ■ 

Why is such result of physical interest? (i) Having con- 
structed the action S gives a closed form solution for the 
path integral, (ii) The action 5' defines a renormalized ac- 
tion in Q.M. (iii) Gutzwiller's trace formula has been 
used widely to study quantum chaos in a semi-classical 
regime (e.g. highly excited states in atoms). Eq.(^) es- 
tablishes a simple relation between the transition am- 
plitude and some action S which mathematically has the 
form of the classical action. Eq.(y) allows to define quan- 
tum chaos in an unambigous way. (iv) Eq.(^ allows to 
study quantum descendents of classical instanton solu- 
tions. We have no proof of this conjecture. However, we 
will give numerical evidence in support of the conjecture. 
We will elaborate on the topic (ii), while (iii,iv) will be 
discussed in another letter. 



where the Euclidean quantum action is given by Se — 
J dtYi^ + y{x)- The transition to imaginary time causes 
a sign change of the mass term, besides that the action 
parameters keep its absolute and relative values. In the 
sequel of the paper we will work in imaginary time. We 
drop the subscript E in what follows. 

We consider a particle of mass m moving in the pres- 
ence of a local potential V{x), given in polynomial form 
Vix) = X)n=o ""2;", with the parameters m,vo, . . . ,vn- 
We write the potential of the quantum action V{x) — 
^^^Q-Cnx", being parametrized by rh,vo, . ■ . ,vn- We 
choose a value for the time T and a set of bound- 
ary points {xi, . . . ,xj}. For all pairs of Xi„, Xfi from 
that set we compute the Euclidean transition amplitude 
Gij = G{xi,Xj,T). We denote the Euclidean action S 
along its classical path Xd by T,ij = ^[xczlli^'.'o ■ ^^~ 
cause Z does not depend on Xin, xji, we subsume it into 



^sub 



t,j - InZ. Then Eq.(|) takes the form 



2. Construction of renormalized action 

We know the quantum action for the harmonic oscillator. 
What does perturbation theory predict for the quantum 
action in the presence of a small anharmonic perturba- 
tion S[x] = /ff f ±2 + !^x^ + XV{x), A « 1? Using 
the saddle point method one obtains up to order 0{\^) 



S = S 



Si[x] = — 
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where G'(i, t') denotes the Euclidean Green's function. 
This shows that the quantum action contains terms be- 
yond that occuring in the classical action. 

In the following we want to determine numerically the 
renormalized action, Eq.(y), for the case of a quartic po- 
tential. In this work, we search for renormalized parame- 
ters, neglecting potential terms higher than fourth order 
in a first step. This search requires to calculate the tran- 
sition amplitude, Eq.(|]). We have chosen to compute the 
path integral via Monte Carlo with importance sampling. 
This requires to go over to imaginary time (Euclidean 
path integral) t -^ —it. Then Eq.(||) becomes 



GE{xfi,tfi;Xr 



[dx] exp[--5_E[x]] 



Xinjtij^ 



where the Euclidean classical action is given by Se 



JdV. 



V{x). Correspondingly, Eq.(y) becomes 
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G,, = exphS?""], *,j = l 



(5) 



While the quantum action has A^ + 2 parameters, Eq.(0) 
represents J'^ /2 + 0{J) independent equations to deter- 
mine the quantum action. We have chosen a larger num- 
ber of equations than parameters (over-determined) re- 
quiring that the error in Eq.(0) becomes globally mini- 
mal. I.e., we make a x^"fit. 
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which is a function of the parameters of the quantum 
action. A solution for the quantum action is considered 
as consistent, if we find the same action parameters for 
different sets of boundary points. 

3. Numerical results 

We have computed the Euclidean propagator in two 
ways: (i) We use Monte Carlo with importance sampling. 
The trick is to write the propagator matrix element as a 
ratio of two path integrals. This corresponds to splitting 
the action S* = 6*0 -f- 6*1, chosen such that the path inte- 
gral for 5*0 is analytically known, exp[— S'o/?i] is treated 
as weight and exp[— S*!/?!] is treated as observable, (ii) 
According to the Feynman-Kac formula, the propagator 
for large time is asymptotically given by the ground state 
contribution. We have solved the stationary Schrodinger 
equation, computing a number M (of lowest lying) en- 
ergies and wave functions, and expressed the propagator 
as a spectral sum over those states. 

(a) To test our algorithms, we considered the har- 
monic oscillator V{x) = V2x'^ ^ V2 — 1/2 (with param- 
eters m = lli ^ h ^ 1). The classical path, the transition 
amplitude and the renormalized action are analytically 



known, giving S ~ S. The results are found within sta- 
tistical errors to be consistent with S — S. 

(b) Next we have considered the quartic interaction 
V{x) — U4X'*, V4 = 1, with parameters m = h ~ 1. 
We have computed numerically the ground state energy, 
Egr = 0.667986, and the variance of the ground state 
wave function, vavgr ~ 0.287333. From those we intro- 
duce dynamically a time scale and a length scale (cor- 
responding to the Bohr radius), T^c = h/Egr — 1.497, 
Lsc — y/fciTgr — 0.5360. We have computed the prop- 
agator for < T < 2 via Monte Carlo runs with 
Ntherm = 4000, N^Mp = 2000 and N.onf - 1000. For 



T > 1 we also computed the propagator by solving the 
stationary Schrodinger equation and calculating the wave 
function and energy of the lowest M states of the spec- 
trum. We used M = 30 for 1 < T < 2 and M = 7 
for T > 2. We used up to 20 000 mesh points to solve 
the classical equation of motion in the determination of 
5[xc/]. For the time interval T = 0.5, the numerical re- 
sults are shown in Tab.[l]. The following observations 
are made: (i) The linear and the cubic renormalized po- 
tential term are compatible within statistical error with 
the value zero, expected from parity conservation of the 
quantum system. 



Tab.l. Quartic potential V{x) = x'^. Harmonic oscillator weight factor. T = 0.5, J = 6. 





m 


Vo 


Vl 


V2 


«3 


Vi 


x' 




Fit MC 


0.9936(3) 


1.1711(18) 


0.000(7) 


0.449(18) 


0.000(16) 


0.982(23) 


21.35 


[-1.0, +1.0] 


Fit MC 


0.9938(2) 


1.1665(19) 


0.000(7) 


0.488(15) 


0.000(13) 


0.954(15) 


17.47 


[-1.2,+1.2] 


FitMC 


0.9941(2) 


1.1695(20) 


0.000(7) 


0.466(12) 


0.000(11) 


0.975(10) 


24.98 


[-1.4,+1.4] 


Fit MC 


0.9940(2) 


1.1696(21) 


0.000(8) 


0.469(10) 


0.000(10) 


0.973(7) 


63.01 


[-1.6, +1.6] 


Fit MC 


0.9944(2) 


1.1684(23) 


0.000(8) 


0.458(9) 


0.000(9) 


0.984(5) 


38.82 


[-1.8, +1.8] 


Fit MC 


0.9942(2) 


1.1679(25) 


0.000(9) 


0.459(9) 


0.000(9) 


0.986(4) 


55.73 


[-2.0, +2.0] 


FitMC 


0.9939(2) 


1.1721(26) 


0.000(9) 


0.449(8) 


0.000(8) 


0.992(4) 


57.80 


[-2. 2, +2. 2] 


Fit MC 


0.9938(2) 


1.1689(27) 


0.000(10) 


0.460(9) 


0.000(8) 


0.987(3) 


43.45 


[-2.4,+2.4] 


Fit MC 


0.9942(3) 


1.1664(29) 


0.000(11) 


0.447(9) 


0.000(8) 


0.993(3) 


55.34 


[-2. 6, +2.6] 


Fit MC 


0.9938(3) 


1.1738(29) 


0.000(12) 


0.432(10) 


0.000(8) 


0.998(3) 


58.76 


[-2. 8, +2.8] 


Fit MC 


0.9946(3) 


1.1705(31) 


0.000(13) 


0.435(10) 


0.000(9) 


0.990(3) 


44.48 


[-3.0, +3.0] 


Fit MC 


0.9940(4) 


1.1702(34) 


0.000(13) 


0.461(11) 


0.000(9) 


0.984(3) 


64.68 


[-3. 2, +3. 2] 


Fit MC 


0.9943(4) 


1.1645(37) 


0.000(15) 


0.460(11) 


0.000(10) 


0.985(4) 


63.00 


[-3.4,+3.4] 


FitMC 


0.9949(5) 


1.1605(41) 


0.000(16) 


0.483(12) 


0.000(11) 


0.973(4) 


47.03 


[-3. 6, +3.6] 


Average 


0.9941(3) 


1.1685(27) 


0.000(10) 


0.458(11) 


0.000(10) 


0.983(7) 


46.85 





(ii) The renormalized action generates a quadratic po- 
tential term, absent in the classical action. With respect 
to variation of the interval, most coefficients vary very lit- 
tle and are compatible within statistical errors with each 
other. The quadratic term is more sensitive, showing 
fluctuations from its average value of up to two standard 
deviations. In general we find for J = 6 that the numer- 
ical values of x^ lie in the order of the value x^ = 36, 
which means that the fits are acceptable. We found that 
the intervals, where the algorithms work, lie in the range 
of [—1.0, +1.0] to [—3.6, +3.6]. If the interval is too small, 
with boundary points Xi << 1, then the potential terms 
x^ and x** are difhcult to distinguish. On the other hand, 
if the interval is too large, with boundary points Xi » 1, 
then the observable exp[— x**] (of the Monte Carlo algo- 
rithm) becomes very small. In the renormalized action 
the quartic term dominates strongly over the quadratic 
term, which makes the latter difficult to discriminate. 
Stability of the renormalized parameters under variation 
of the interval of boundary points has been observed also 
for other values of the time T. The dependence on T is 
shown in Fig. [1] . The error bars are estimated systemat- 
ical errors (statistical errors of Monte Carlo data are of 



the size of symbols). We distinguish three regimes: (i) 
< T « T,„ (ii) T « T,„ and (iii) T » T,,. For 
T close to zero, the renormalized parameters are close 
to those of the classical action. An exception is V4 at 
T = 0.1, which is about 10% off the classical value. We 
believe that this does not reflect physics but is a fault due 
to limited numerical precision (recall that the propagator 
becomes S{xfi — Xm) when T —^ 0). At about T ss Tsc, 
rh, V2 and W4 change noticeably. For large T >> Tsc the 
renormalized parameters converge asymptotically. When 
increasing T beyond Tgc one observes that the stability of 
results requires an exponential increase in the number of 
mesh points. This can be traced to the fact that the clas- 
sical trajectory covers several orders of magnitude (and 
the order increases with T). This observation, similar to 
critical slowing down observed in simulating critical phe- 
nomena in lattice field theory, puts an upper limit on the 
time parameter used in this investigation (T < 5). 



4. Interpretation 

How can we understand the behavior of the renormalized 
parameters? Why does the renormalized action depend 
on T? Consider the n-point function of a scalar field. 



r =< vac\t.o.{^{xi) ■ ■ ■ ^{xn)\vac > . 

This is the vacuum expectation value of a time-ordered 
n-fold product of the field. It corresponds to the transi- 
tion from the physical vacuum ait — — cxd, to the physical 
vacuum at i = -|-oo with intermediate creation (annihi- 
lation) of particles at ii, • • • , t„. In order to obtain finite 
expressions, one introduces a regularisation parameter, 
say /i. Then F = T{xi,fj.), i.e. the result depends on 
/i. In Q.M., we consider the propagator from Xi„,ii„ to 
Xfi, tfi, T ^ tfi — tin- The analogue in Q.F.T. is 

< vac,t = r|t.o.($(xi) • • ■^{xn)\vac,t = > . 

Then the n-point function T(xi, fi.T) depends also on 
the time T. The computation of renormalized parame- 
ters from T{xi^ /z, T) then gives mass, coupling constants 
etc. as function of T. Thus it is no surprise that the 
renormalized parameters in Q.M. also depend on T. 



G{xfi,T;Xin,0) '^T^oc ll^gr{xfi) exp{-EgrT/h)'lljgr{xin). 

For the harmonic oscillator {Egr = huj/2) one can show 
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n^OO E, 



gr- 
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For the quartic potential one observes also a smooth be- 
havior of vq. We have fitted vq by the function A + B/T. 
We obtain vq -^ 0.6686, compatible with Eq.(^. Then 
using the definition of the renormalized action, Eq.(0), 
splitting the renormalized action into a vq part and a 
rest, and using the Feynman-Kac formula implies 



>T^oo ipgrixfi) e:iq){-EgrT/h)ipgr{x^, 



V{Xcl)] 



Consequently, the renormalized action, when excluding 
the Vq term, has an asymptotic limit for large T 



1 I ^ Tfi • 2 



h\[i:gr{xfi) IpgriXin)]- 



This is a strong indication (no proof) that the renormal- 
ized parameters of S converge asymptotically in T . This 
is what has been observed in the data. It corresponds in 
Q.F.T. to the infinite volume limit. 

In conclusion, we have computed numerically for the 
quartic potential a quantum corrected action, and find 
parameters being quite different from the classical val- 
ues. Hence, quantum fluctuations can change the clas- 
sical behavior drastically. This has been observed (in 
preliminary results) also for the instantons solutions of 
the double well potential. 



Fig.[l] Renormalized parameters (in dimensionless units) 
versus time T . Transition amplitude obtained by Monte 
Carlo (dots and circles) and by stationary Schrodinger 
equation (other symbols). 
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How can we understand the behavior of the renor- 
malized parameters changing qualitatively in the inves- 
tigated time interval? (i) For sufficiently small T, one 
has S < h, i.e. we are in the quantum regime. However, 
according to Dirac ||§] (see also Rcf. jlj), ex.p{iS/h) is a 
good approximation of the propagator G, when the time 
interval T over which G is supposed to propagate goes to 
zero. This is consistent with our numerical observation 
S ~ S in this regime, (ii) For sufficiently large time T 
one has S > h, i.e. one is in the semi-classical (WKB) 
regime. In imaginary time, the Feynman-Kac formula 
describes the propagator asymptotically. 



[1] L.S. Schulman, Techniques and Applications of Path In- 
tegration, John Wiley&Sons, New York (1981). 

[2] L.S. Schulman, Phys. Rev. A35(1987)4956. 

[3] L.S. Schulman, Phys. Rev. 176(1968)1558. 

[4] Private communication by E. Nelson. See |l[. 

[5] M.C. Gutzwiller, Chaos in Classical and Quantum Me- 
chanics, Springer, Berlin (1990), and references therein. 

[6] P.M. A. Dirac, Physikalische Zeitschrift der Sowjetunion, 
3, No. 1 (1933); reprinted in: J. Schwinger, Quantum Elec- 
trodynamics, Dover, New York (1958); see also P.M. A. , 
The Principles of Quantum Mechanics, Oxford, London 
(1958). 



1.2 



1.0 - • i 



0.8- 



0.6- 



0.4- 



▲ 

A 



o 


M 


.C. 


ik 


S. 


E. 30 states 


A 


S. 


E. 7 states 




0.2- 





1.0 H o o o o 

0.8 
1.0-1 



it 



m 



V, 



A 



A 



0.8- 



A 



A A 



■gi" 



0.6 



"SC 



2 



T 
3 



~r 

4 



T 
5 



